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Abstract 

The production of a boson and two jets, at least one of which contains a b quark, 
is a principal background to single-top production, Higgs production, and signals of 
new physics at hadron colliders. We present a next-to-leading-order (NLO) calculation 
of the cross section at the Fermilab Tevatron and the CERN Large Hadron Collider. 
The NLO cross section differs substantially from that at LO, and we provide a context 
in which to understand this result. 



1 Introduction 



Many signals for new physics at hadron colliders involve an electroweak gauge boson (7, Z, W) 
and jets, one or more of which contain a heavy quark (c, b). The prime example is the dis- 
covery of the top quark via the signal W + Aj, where two of the jets contain b quarks [H [2]- 
It is important to understand the backgrounds to these signals in as much detail as possible 
[3]. We have recently completed a next-to-leading-order (NLO) calculation of Z production 
in association with one or two jets, one (or more) of which contain a heavy quark |5J, |6j 111 
In this paper we present a NLO calculation of W production in association with two jets, 
one or more of which contain a b quark. The case of W production in association with one or 
more jets containing a charm quark is more complicated due to the additional contributions 
coming from s — > Wc. 

The production of a boson in association with two jets, one or more of which contain 
a heavy quark, is particularly interesting as it is the principal background to both single-top 
production [H [9l [10] and Higgs production (via qq' — > Wh) [HI [12], which are currently 
being sought at the Fermilab Tevatron. Single-top production occurs via both a t-channel 
process, qb — > q't [TSj [HJ [15], and an s-channel process, qq' tb [HI [T7]. In both cases 
the final state, after top decay, comprises a W boson and two jets, one (t-channel) or both 
(s-channel) of which contain a b quark. In the case of Higgs production, the final state (after 
the decay h bb) comprises a W and two jets, both of which contain a b quark [THl 
However, even if both jets contain a b quark, it is more efficient to tag only one of them 
|12j . Thus in all cases the principal background is W plus two jets, one (or more) of which 
contain a heavy quark. The signals for single top[20l[ni[2a[2l[2ll25llMl[27l[2Hl[29l[30] 
and Higgs [311 1321 EH] are known at NLO and beyond; our goal is to provide a calculation of 
the principal background at NLO. 

At leading order (LO), there are two processes that produce a W boson and two jets, 
at least one of which contains a b quark. These are shown in Fig. [H The process qq' — > 
Wbb yields two b jets, and is already known at NLO [Ml [351 [26]. In contrast, the process 
bq — >■ Wbq' yields just one b jet. For a signal with just one b tag, these two processes are 
comparable at the Tevatron, while at the CERN Large Hadron Collider (LHC) the latter 
process is dominant [37] . 

The process bq Wbq' requires further consideration, as it contains a b quark in the 
initial state. The b distribution in the proton is derived perturbatively via the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations [5S1 [SH]- Alternatively, one 
may eschew a b distribution function, and use a gluon in the initial state, which then splits 
to bb, with one b participating in the hard scattering while the other remains at low transverse 
momentum; this is shown in Fig. [2] [37J . However, there are two advantages to working with 
a b distribution function. First, initial-state collinear logarithms, of order [ag ln(Miy/mfe)]"', 
are summed to all orders, yielding a more convergent perturbative expansion. Second, the 
LO process is simpler, which makes the NLO calculation tractable. 

A b distribution function is useful when the scale of the process is much greater than the 
b mass. In beauty production at HERA, the scale is not much greater than the b mass, and 
the appropriate description is 7* (7 — > bb [IQ]. Thus there is no direct measurement of the b 

^The inclusive production oi a. Z with one or more heavy quarks is presented in Ref. [7]. 
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Figure 1: Leading-order processes for the production of a boson and two jets, at least 
one of which contains a b quark. 
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Figure 2: An alternative way of calculating the process in Fig. [T](b). 

distribution function at this time. Single-top production (t-channel) may provide the first 
such measurement. 

The paper is organized as follows. In Section [2] we outline the NLO calculation and 
discuss some of the finer points. The issue of correctly treating the 6-quark mass arises in a 
novel way, and we discuss this in some detail. The results of the calculation are presented 
in Section [31 



2 Wbj at NLO 

The LO cross sections for Wbj and Wbb, arising from the processes in Figs. [11(b) and[Tl(a), 
respectively, are given in parentheses in Table [H All jets are required to have px > 15 GeV 
and \r]\ < 2 at the Tevatron, while at the LHC we require px > 25 GeV and \ri\ < 2.5. The 
jets are required to be separated by ARjj > 0.7. As mentioned in the Introduction, these 
cross sections are comparable at the Tevatron, while Wbj is dominant at the LHC. This 
qualitative result is already well known from the LO calculation of Ref. [37], which used a 
gluon in the initial state (Fig. [2|) rather than a b quark [Fig. [H^b)]. Our goal is to calculate 
these cross sections at NLO. 

The processes involved in the calculation are as follows: 
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qq' - 


Wbb at tree level [Fig. 


[T](a)] and one loop 


• 


bq - 


^ Whq' at tree level [Fig. 


[T](b)] and one loop 




nn^ - 
HH 
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• 


bq - 


^ Wbq'g at tree level 




• 


m - 


->■ Wbbq' at tree level (Fii 




• 


bg- 


Wbqq' at tree level 





In all cases q = u,d, s, c, and we include Cabbibo-Kobayaslii-Maskawa (CKM) mixing. The 
6-quark mass is neglected throughout (except where noted). The calculation is simpler than 
the corresponding calculation for ZQj and ZQQ because there are several subprocesses that 
have no contributing analogue, namely Qg ZQg and gg ZQQ [6J. The analogous 
processes, bg Wtg and gg Wtb, yield a top quark in the final state. We do not consider 
these, or any other processes that yield a top quark, as part of our calculation of the Wbj 
and Wbb backgrounds. 

The NLO calculations in this paper were performed with the Monte-Carlo code MCFM 
H2] • The leading-order calculations were performed both with this code and with MadE- 
vent [l3]. The NLO corrections are included in the code MCFM by implementing the virtual 
helicity amplitudes of Ref. [H]. These are given in the four-dimensional helicity scheme, 
which is used throughout the calculation. The strong coupling constant is shifted from this 
scheme to the MS scheme The real corrections are adapted from Ref. [16] and sin- 

gularities are handled using the dipole subtraction method [17] . Since the matrix elements 
contain the decay of the gauge boson into massless particles, the code is general enough to 
provide results for final states from W* — *• iu, including the correlation between the W* spin 
and the angular distribution of the leptons. For this paper we specialize to the case of an 
on-shell W boson. 

The NLO cross sections are presented in Table [H To obtain the NLO cross sections for 
Wbj and Wbb, which involve the radiation of additional partons, two partons are combined 
into a single jet by adding their four-momenta if ARjj < 0.7. If the combined partons are 
both b quarks, then the process contributes to the column labeled W{bb)j. This is a.W + 2j 
event in which one jet contains two b quarks, which changes the tagging probability for that 
jet [IE]. It is calculated with a finite 6-quark mass {rrih = 4.75 GeV) in order to regulate the 
logarithmic divergence present when the heavy quarks are coUinear. If all three partons are 
well separated, then the process contributes to either Wbbj or Wbjj. 

There is another aspect of the calculation in which the 6-quark mass cannot be neglected. 
If the b quark comes from a final-state virtual gluon splitting to bb, with the other b quark 
missed, then the cross section is sensitive to the 6-quark mass. This can occur in the NLO 
processes qq' —>■ Wbbg and gq —>■ Wbbq', both of which contribute to Wbj when one b quark 
is missed. We therefore calculated these processes with a finite 6-quark mass. In the case of 
the process gq — >■ Wbbq', which involves the splitting g ^ bb in the initial state (see Fig. [2]) 
as well as the final state, this requires that we abandon the dipole subtraction method in 
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favor of subtracting the mass singularity via the truncated b distribution function [381 ES] 

where Pqg{z) = ^[z"^ + (1 — 2)^] is the DGLAP sphtting function. The counterterm is 
constructed by calculating bq Wbq'; this cancels the initial-state logarithmic dependence 
on the 6-quark mass in gq Wbbq', and yields a cross section in the MS factorization 
scheme. 

It is interesting to compare the LO result for Wbj with a calculation based on gq Wq'bb 
(Fig. [2]) [37|, which does not use a b distribution function. With the same parameters as 
Table m and using rrij, = 4.75 GeV and CTEQ6L1 parton distribution functions [49], we 
find a cross section of 1.92 pb at the Tevatron. This is to be compared with 1.06 pb using 
bq — *■ Wbq', and serves as a rough check of the formalism. Both of these numbers are LO 
and very scale dependent, so one should expect only qualitative agreement H The result at 
the LHC {W~^ + W~) is 81 pb, which is to be compared with the 87.3 pb result obtained 
using bq — > Wbq' (see Tabled]). In any case, the NLO results presented in this paper should 
be more accurate than any of these LO results. 

We checked that the effect of the heavy-quark mass is small if all b quarks are at high pt 
by comparing Wbb at tree level with and without a finite quark mass. With the cuts used 
in this paper, we find a reduction of the Wbb cross section by about 7% at the Tevatron 
{Pt > 15 GeV) and 2% at the LHC {pr > 25 GeV). This is the same as the result found 
for the effect of the 6-quark mass on the full NLO calculation at the Tevatron [3^], which 
suggests that the effect of the 6-quark mass can be divined from a LO calculation. We found 
that the heavy-quark mass is a similarly small effect for Wbbj at tree level, as expected. 

We also list in Tabled] the LO (in parentheses) and NLO cross sections for Wjj [35] [50]. 
and the LO cross section for Wjjj. We ignore CKM mixing, which is a negligible effect 
at LO, and is therefore also negligible at NLO. In these cross sections we have included 
the contribution from light partons as well as heavy quarks (c, 6). Thus, for example, the 
fraction oi W + 2j events in which only one of the jets can be tagged as b jet is given by 
[Wbj + W{bb)j]/Wjj. 

3 Results 

The results in Tableware remarkable. First we see that, at NLO, the cross sections for Wbj 
and Wbb are comparable at the Tevatron. Thus both final states are important backgrounds 
to single top and Higgs. At the LHC, Wbj is much larger than Wbb. At both machines, the 
correction to Wbj is large (for fip = f^R = M^), about a factor of 1.9 at the LHC and a 
factor of 2.4 at the Tevatron. The correction is larger at the Tevatron because Wbb, which is 
comparable in size, feeds into the Wbj column at NLO when one of the b quarks is outside 
the fiducial region {pt > 15 GeV, \ri\ < 2). 

^Including the uncertainties due to the independent variation of the renormahzation and factorization 
scales from Mw/2 to 2Mw, the cross sections are l-92tQit-QA3 pb for gq Wq'bb and l-OetoJo to.o7 Pb 
for bq Wbq' 
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The correction to Wbb is significant at both machines, but more modest than that of 
Wbj, again for = fi^ = M^. At the LHC the cross section for Wbbj is comparable to 
that of Wbb, which calls into question the reliability of perturbation theory. However, Wbbj 
is also a correction to Wbj, and when compared with that LO cross section it is a modest 
correction. Nevertheless, it does mean that when demanding two b quarks in the final state, 
the processes Wbb and Wbbj are comparable at the LHC. The latter is only known at LO and 
therefore has a large uncertainty. This is reflected by the large scale dependence associated 
with this process, as discussed in Refs. [35| [36]. 

Using the results in Table [H we construct a set of inclusive cross sections, presented in 
Table O Displayed this way, the radiative corrections to the cross sections are even more 
significant, especially for Wbb + X at the LHC, where the NLO cross section is about 2.7 
times the LO cross section. However, the proper way to regard this is that a new process 
enters at NLO, namely gq Wbbq' (see Fig. [2]), which is also a correction to the much larger 
LO process bq — > Wbq' , as discussed in the previous paragraph. We also list the uncertainties 
from varying the renormalization scale from My// 2 to 2M\jy while keeping the factorization 
scale fixed at My/, from varying the factorization scale from My// 2 to 2My/ while keeping 
the renormalization scale fixed at My/, and from varying the parton distribution functions 
|49j . The largest uncertainty is from varying the renormalization scale, since the process 
under consideration is of order q;| at leading order. 

We show in Figs. [3] and H] the renormalization- and factorization-scale dependence at LO 
and NLO for the inclusive Wbj +X cross section at the Tevatron and the LHC, respectively. 
At the Tevatron, there is almost no reduction of the renormalization-scale dependence at 
NLO. This supports our earlier argument that part of the source of the large NLO correction 
to this process is that the process Wbbj contributes when one of the b jets is missed. Since 
this is a tree-level process, it has a large renormalization-scale dependence. There is only a 
mild reduction of the renormalization-scale dependence at the LHC. At both machines there 
is a mild reduction of the factorization-scale dependence. 

Using the code MCFM [HI |12], we are able to plot any desired distribution at both LO 
and NLO. For example, we show in Fig. [5] the exclusive NLO differential cross section at 
the Tevatron for Wbb, Wbj, and W{hb)j, versus the transverse momentum {pt) of the jet 
with the highest pt- For Wbj, the highest-p^ jet contains a b quark 48% of the time; for 
W{bb)j, the percentage is 63%. The W{bb)j cross section, which first arises at NLO, has a 
very different shape from the other two, and is only significant at high px- 

The lower histogram in Fig. [5] shows the ratio of the NLO and LO cross sections. We see 
that there is a change in shape of the Wbj and Wbb differential cross sections at NLO, with 
relatively more events at low px- There is also an enhancement of the Wbj cross section at 
high Pt, due to the process qq — > Wbbj with one b missed. 

We show in Fig. [6] the exclusive NLO transverse momentum distribution of the charged 
lepton from W decay at the Tevatron. We again see a change in shape at NLO, with relatively 
more events at low pt for Wbb, and relatively more events at high p^ for Wbj, the latter 
again due to the process qq — > Wbbj with one b missed. This plot demonstrates the abihty of 
MCFM to calculate quantities involving the decay products of the W boson while correctly 
treating the correlation between the angular distribution of the decay products and the W 
spin. The W boson is treated as on-shell in these plots, but MCFM is flexible enough to 
allow for off-shell W bosons as well. 
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Table 1: Exclusive cross sections (pb) for W boson plus two (or more) jets, with at least 
one b jet, at the Tevatron {^/s = 1.96 TeV pp, pr > 15 GeV and \r]\ < 2) and LHC 
{^/s = 14 TeV pp, pt > 25 GeV and |?7| < 2.5). Two final-state partons are merged into 
a single jet if ^Rjj < 0.7. No branching ratios or tagging efficiencies are included. The 
labels on the columns have the following meaning: Wbj = exactly two jets, one of which 
contains a b quark; Wbb = exactly two jets, both of which contain a b quark; W{bb)j = 
exactly two jets, one of which contains two b quarks; Wbjj = exactly three jets, one of 
which contains a b quark; Wbbj = exactly three jets, two of which contain a b quark. For 
the last set of processes, which include both light and heavy partons in the final state, the 
labels mean: Wjj = exactly two jets; Wjjj = exactly three jets. For Wbj, Wbb, and Wjj, 
both the leading-order (in parentheses) and next-to-leading-order cross sections are given. 
The CTEQ6M parton distribution functions are used throughout, except for the LO cross 
sections in parentheses, where CTEQ6L1 is used [49j. The factorization and renormalization 
scales are chosen as fip = fJ^R = M^. 





Cross sections (pb) 


Collider 


Wbj Wbb W{bb)j 


Wbjj Wbbj 


TeV W+{= W-) 


(1.06) 2.54 (2.48) 3.14 0.89 


0.18 0.65 


LHC W+ 


(51.7) 96.2 (9.5) 14.3 27.0 


13.8 11.6 


LHC W^ 


(35.6) 66.4 (6.6) 9.6 19.0 


9.3 7.6 




Wjj 


Wjjj 


TeV W+{= W-) 


(261) 290 


39 


LHC W+ 


(4990) 4170 


1280 


LHC W' 


(3650) 3030 


890 



We show in Fig. [7] the exclusive NLO di-jet invariant mass at the Tevatron. Both Wbb 
and Wbj are enhanced at low invariant masses at NLO. The process W{bb)j is significant 
at high invariant mass. All three processes are numerically important as backgrounds to the 
search for a Higgs boson via the process qq Wh — Wbb. For example, the integrated 
cross sections in the region 110 GeV < rrijj < 130 GeV are 130 fb for Wbj, 95 fb for Wbb, 
and 75 fb for W{bb)j. 

We show in Figs. lHM fTOlthe same distributions at the LHC, but this time for the inclusive 
cross sections, rather than exclusive. At the LHC, the process Wbj -f- X is dominant, while 
Wbb + X and W{bb)j are comparable. There is little change in shape of the distributions 
at NLO, with the exception of the px spectrum of the highest pt jet, which is enhanced at 
high Pt for Wbb + X. The di-jet invariant mass distribution also shows an enhancement at 
low values of the invariant mass, but not at values relevant for the Higgs search. 

The code MCFM is publicly available [12], and one can generate any desired distribution. 



6 



Table 2: Inclusive cross sections (pb) for W boson plus two (or more) jets, with at least one 
b jet, at the Tevatron {^/s = 1.96 TeV pp, px > 15 GeV and |?7| < 2) and LHC {^/s = 14 
TeV pp, pt > 25 GeV and |?7| < 2.5). Two final-state partons are merged into a single jet 
if ARjj < 0.7. No branching ratios or tagging efficiencies are included. The labels on the 
columns have the following meaning: Wbj + X = at least two jets, one of which contains a 
b quark; Wbb + X = at least two jets, two of which contain a b quark; W{bb)j = two jets, 
one of which contains two b quarks. The label Wjj + X = at least two jets, containing both 
light and heavy partons. For Wbj + X, Wbb + X, and Wjj + X, both the leading-order 
(in parentheses) and next-to-leading-order cross sections are given. The CTEQ6M parton 
distribution functions are used throughout, except for the LO cross sections in parentheses, 
where CTEQ6L1 is used [19]. The factorization and renormalization scales are chosen as 
A^F = A^R = Mw- The uncertainties are from the variation of the renormalization scale, the 
factorization scale, and the parton distribution functions, respectively. 





Cross sections (pb) 


Collider 


Wbj + X Wbb + X 


Wibb)j 


TeV W+{= W-) 
LHC W+ 
LHC W- 


n riR^ 79+0.68 +0.01 +0.09 (n aQ\ Q 7Q+0-56 +0.12 +0.11 
U-UDJ A. ( Z_o.50 -0.04 -0.09 l^-^yj O. (i)-0A7 -0.08 -0.11 

(51.7)1101111^^11 (9.5) 25.911? 1° 3 lg| 
(35.6) 75.7lL^l^:il|^ (6.6) 17.2111 l°Jl°:i 


n SO+O-36 +0.10 +0.03 
U-oy-0.23 -0.07 -0.03 

97 n+11 +0-1 +0-9 
^ ' -U-T.l -0.3 -0.9 

1Q n+^-^+o-i +0-7 

iJ.U_5,o -0.3 -0.7 




Wjj+X 


TeV W+{= W') 
LHC W+ 
LHC W- 


(261) 32911° 1^1? 
(4990) 5450111^11' lll'o 
(3650) 3920ll?r?ritS 



We hope the calculation outlined in this paper is useful to better understand the backgrounds 
to single-top production, Higgs production, and signals of new physics at the Tevatron and 
the LHC. 
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Figure 3: Renormalization- and factorization-scale dependence of the inclusive Wbj+X cross 
section at LO (dashed) and NLO (solid) at the Tevatron. The cross section is normalized to 
its value at the reference scale /xq = Mw- 
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Figure 10: Same as Fig. [HI but for the di-jet invariant mass. For the invariant 

mass is that of the two b jets; for the invariant mass is that of the b jet and the 
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15 



